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' Abstract 

' The SchifF moment of the ^®^Hg nucleus was calculated using finite 

range P- and T-violating weak nucleon-nucleon interaction. Effects of the 
core polarization were considered in the framework of RPA with effective 
residual forces. 

> 

oo 

^ ' 1 Introduction 

' The most precise limit on parity- and time-invariance violating nucleon-nucleon 

, interaction has been obtained from the measurement of the atomic electric dipole 

' moment of ^^^Hg 1 . The hadronic part of the atomic dipole moment associated 

with the electric dipole moment of the ^^^Hg nucleus manifests itself through 
the Schiff moment which is the first nonzero term in the expansion of the nuclear 
electromagnetic potential after including the screening of the atomic electrons 

[213111. 



I The operator for the Schiff moment is [H] 



where is e for a proton and zero for a neutron. The Schiff moment generates 
T- and P-odd electrostatic potential in the form 

0(r) = 47rS • VS{r). (2) 

Interaction of atomic electrons with the potential given by Eq. |2l produces an 
atomic dipole moment 



-eEf<l>ir^)\n){n\~eJ:fzM ^ ^^^^ 

En — Eq 
ft. 

Due to contact origin of the potential, the electrons in s- and p- atomic orbitals 
only contribute to the dipole moment given by Eq. Q- 
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The nuclear SchifF moment was calculated so far in a simplified model plE] 
without considering many-body nuclear structure effects. These effects have 
to be understood properly if we intend to extract the parameters of P- and 
T-violating nuclear interaction from the value of the Schiff moment. For light 
nuclei the properties of the Schiff moment strength obtained in modern shell 
model calculations were discussed in Ref. |7j. The study of the polarization 
effects associated with the coupling to isoscalar dipole compression mode has 
been performed in Ref. [5]. In our paper we calculate the Schiff moment of 
is^Hg nucleus within RPA framework with effective residual strong forces using 
finite range P- and T-odd weak nuclear interaction. 



2 Basic ingredients of the theory 

2.1 Nucleon-nucleon P- and T-odd interaction 

We use the interaction generated by P- and T-violating pion exchange |9I10I1H . 

W{ri - r2) = - ^ [{goTi ■ T2 + 52(ti ■ T2 - 3rf r|))((Ti - 0-2) 

+9i{t^o-i - tI(T2)] ■ Vi , (4) 

ri2 

where g is the usual strong pion-nucleon pseudoscalar coupling constant, go, gi, 
and g2 correspond to isoscalar, isovector , and isotensor P- and T-odd couplings, 
nip is the proton mass. In contrast to P-odd and T-even interaction, in Eq. Q 
the exchange of tt" is allowed. This term produces the direct contribution to 
P- and T- odd part of nuclear mean field while the other terms produce the 
exchange contribution only. Since the direct contribution dominates for finite 
range potentials we can expect that the interaction 10} is the leading one and 
the exchange of heavier mesons can be omitted. 

In previous calculations the phenomenological contact interaction has often 
been used instead of finite range interaction given by Eq. 10}. It has the form ^2] 

G 1 

Wcira - Vb) = "^^^^ iiVabfTa - VbatTb) ' ^a^iVa - rb) + 

Vabif^a X CTb] • {(p^ - Pfc), (5(ra - Tfe)}) , (5) 

where G is the Fermi constant. In the limit m^r —> oo the interaction (0} 
transforms into Eq. (|S} after the substitution ggi — s- ^ rj. We shall use 
this factor when comparing our results with those obtained using the contact 
interaction given by Eq. 



2.2 Nuclear mean field and correction from the weak forces 

In our calculations we used full single-particle spectrum including continuum. 
The single-particle basis was obtained using partially self-consistent mean-field 
potential of Ref. [El- The potential includes four terms. The isoscalar term is 
the standard Woods-Saxon potential 

U^{r) = - ^ (6) 
1 + exp ^ 
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with the parameters V = 52.03 MeV, R = 1.27Q%A^/^ fm, and a = 0.742 fm. 
Two other terms Uis{r), and Ur{r) were calculated self-consistently usmg two- 
body Migdal-type interaction of Ref. ^1] for the spin-orbit and isovector parts 
of the potential. The last term is the Coulomb potential calculated for uniformly 
charged sphere with Rc = 1.18^-'^/'^ fm. The mean field potential obtained in 
this way produces good fit for single particle energies and r.m.s. radii for nuclei 
in the lead region. 

The correction to the mean field © from the weak interaction Q consist of 
direct and exchange terms. The direct term has the form 

SUdUr) = ^{a- ■ ny / r'^dr' 6io(r,r')[(.go - 2g2)(PpW - p„(r))+ 
Trmp Jo 

giiPpir) + Pn{r))l (7) 

where the function feio(?'i,''2) is a combination of spherical Bessel functions of 
imaginary argument 

hii2{'ri,r2) = ih{m^ri)ki^{m^r2)9{r2-ri)~ii.AmT,r2)kiAm^ri)e{ri-r2). (8) 

Note, that the potential given by Eq. O is pure isovector. The contribution 
of the isovector interaction component dominates in Eq. |(7J, the isoscalar and 
isotensor components of the interaction are suppressed by the factor ■ For 
zero range interaction the potential would be proportional to Vp(r). The gra- 
dient makes this potential very sensitive to the details of nuclear surface. Our 
potential given by Eq. iQ is less sensitive to the surface due to additional 
integration over region of the order of pion Compton wavelength. 

The exchange term is more complicated. The matrix element of it taken 
over angular variable is a non-local operator in radial coordinates. 

{v\5Ue^,h{r,v')\v) = W,{ry), 

where \v) = — (cr • n) 1 1/) , and 



2j. 

h h 1 




bi,i,{r,r')R^{r)R^{r') 



(-)'ny(i>||rUI'^)H|r/,1|«:)*]}. (9) 

The Trace is assumed over isospin variable of the second particle , is the 
occupation number of the single particle state [I] = ^/2l + 1 , and the 
tensor operator Tjj^(n) — {cr ® YL{n)}jM- 

The correction to the single particle wave function ■01/ (r) can be presented 

as 

5V.(r) = {(T-n)^,{n)SR,{r), 

where is the angular part of the wave function. The radial correction 

5Ri^{r) is the sum of the direct and the exchange terms 5R^{r) = SR^dirif) + 
5Ruexch{r), where 

/•oo 

SRud^r{r)^ / G,lAr,r'\e,)5Ud^r{r')RAr'y^dr', 
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SR,e.ch{r)^ G^jy,r'\e,)W,{r' ,r")R,{r"y\"^ dr'dr" . (10) 

Here G- i {r,r'\e,j) is the Green function of radial Schrodinger equation for the 
total angular momentum j^, and the orbital angular momentum = — l^, . 
is the single particle energy. 

3 Core polarization 

The effects of the core polarization for a one particle operator can be treated 
introducing a renormalized operator S satisfying the equation 

S.V = S°,, + ^ S^^,^^^i^—^{v'tx'\F + (11) 

where is the bare Schiff moment operator given by Eq. and e^^ are 

the single particle occupation numbers and energies. For static moments the 
external frequency — > 0. The interaction in Eq. includes both the strong 
residual interaction F and the weak one. The latter we take in the form given 
by Eq. Q). Strictly speaking, the interaction of two nucleons in nuclear matter 
differs from the interaction in the vacuum. We do not discuss this effect here and 
keep the weak interaction in the form Q . The single particle wave functions in 
Eq. are the eigenstates of the mean field which is also the sum of the strong 
and the weak fields given by Eqs. H7I9|1 . Since the weak forces are really small 
compared to the strong interaction, it is natural to treat them perturbatively. 
The simplest way to do it is to present <y\v >= V'iy(r) + Sip„{r) and to gather 
the terms linear in ^^p and VF(ri — r2). It is convenient to present the matrix 
element of the Schiff moment as a sum of three terms 

S,, = W.|S|^,) + (V.|S|^^.) + (V.I^SI^,). (12) 

The operator S in the first two terms satisfies the same Eq. (|11|) where the strong 
interaction only is included. The corrections from the weak forces for these terms 
are in the wave functions of an odd nucleon only. The third term presents an 
induced contribution |15j arising from the P- and T-violating corrections to the 
intermediate states and The equation for SS is 

{6S - 6Snl).'. = (<5So).v + ^(^S - JS^vl)^;.' '"'^ ~ W'^i'\F\H, (13) 

where SSnl is the non-local part of <5S produced by the exchange matrix ele- 
ments of weak interaction W . 

(JSnl).'. - ^(S)^^. "^'""^' {,y'ti'\We.ch\H . (14) 

The equation for SSq is 
(<5So).v = V ~ ((/x|S|M')(z.V|Wd„|H + {f,\SSNL\^^'){l^'fI'\F\^„y) + 

{6^^\s\fi'){iy'fi'\F\f„.) + {^i\s\s^p^,){,y'^l'\F\^i,y)+ 
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{^i\s\^Ji'){v'5i,^,\F\^iv) + ti'\F\5i;^v)) . (15) 

Note, that in the absence of the core polarization the first term only in Eq. Ijf 5|l 
contributes into JSq . And this is the only term that produces the Schiff moment 
of the nucleus with an odd neutron, like ^^^Hg. The residual interaction F has 
the form 

F^C (fir) + f'{T, ■ T2) + gs{<Ti ■ <T2) + g'AcTi ■ <T2)(ti • rs)) ^(ri-ra), (16) 

where C =300 MeV/m^, and f{r) = fex + {f^n- fex)^- Note that the angular 
dependence of the operators S and (5S is completely different. While ~ 
Yi^(n), the induced part (55*^ is a superposition of spin dependent operators 

and {a ® l2(n)}ip. For this reason, different parts of the interaction H16(l 
contribute into renormalization of S and (5S. 

For a spherical nucleus we can separate the angular variables and solve the 
obtained equations in coordinate space. The equations are 

/"OO 

5"(r) = S^{r) + / A''\r, r')S\r') dr' , (17) 
Jo 

where a = p,n and S^ir) is the radial part of the Schiff moment operator Eq. 
(|T|) multiplied by r. The particle-hole propagator Air^r') was calculated by 
means of the Green functions of radial Schrodinger equation. 

A{ry) = ^ •Tr2{(/(r) + /'(xi .X2))^n,|(j7||yi||«;)|Vi?,(r)x 

Kjl 

r'R^{r'){G,i{ry\e^ + uj) + Gji{r, r'|e, - uj)}. (18) 

Similar equations can be written for 5S{r) - the local part of the induced mo- 
ment. They differ from Eq. (|18|l in type of tensor operators and in residual 
interaction. 

3.1 Separation of the spurious component 

The integral equation IjlTfl must have zero eigenmode related to the center of 
mass motion. However, in our case the mean field Eq. ((HJ is not consistent with 
the interaction Eq. H16|l and, in general, we do not have zero energy for the 
center of mass motion. The situation can be improved using some freedom in 
the value of the interaction constant /,„ in Eq. H16(l . We can fix the value by 
the condition loq — Q for the lowest isoscalar dipole mode. Since this procedure 
is numerical the condition uj^ — Q cannot be fulfilled exactly but with finite 
accuracy. For ^^^Hg we found that the set fm = 0.3935, fex = —2.6, and 
/' = 1.07 gives for luq the value = 0.1 KeV which is really small compared 
to the energy of dipole transitions. 

The finite accuracy in determination of the spurious mode brings another 
problem for solutions of the Eq. H17() . The bare operator So{r) becomes non 
orthogonal to the spurious mode transition density. This results in admixture 
of the spurious component to all solutions of the Eq. H17|l . The spurious com- 
ponent should be subtracted since it can change the solution considerably. The 
subtraction can be performed using analytical properties of the solution as a 
function of lo. The solution of the inhomogeneous linear integral equation H17|) 
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as a function of lu has the first order poles at lu = ±Eex, where E^x are the 
excitation energies of the RPA modes. Thus, at small oj the solution can be 
presented as 

5(rM = + bir\u; = 0) + O(^). (19) 

The first term in Eq. H19(l is the spurious component contribution. Let us define 
the following set of integrals in a complex tj-plane along a circle with the radius 
satisfying the conditions t^jQ <C \u}\ <C E^x- 

J„{r)^ (fco^S{r\u)^. (20) 

Using Eq. H19|l we can perform the integration analytically As a result we obtain 

6(r|cj = 0) = J_i(r), 

^o = ^3(r)/Ji(r). (21) 

The numerical integration in Eq. (|20|l has been performed using 32-points Gauss 
formula. The integration radius was |a;| = 0.1 MeV. The obtained field b{r\Lu — 
0) was practically insensitive to the integration radius. The value of ujq was 
independent on r in first 5 or 6 digits. The loss of accuracy was noticeable only 
at \uj\ ~ Wo- 

Fig.l shows the result for renormalized S{r) = b{r\uj — 0) with subtracted 
spurious component. It is shown by full line. The dashed line shows the un- 
renormalized component So(r). The effects of the core polarization are not large. 
We found that for ^"^Bi they change the valence proton contribution by ~ 15% 
which is in fair agreement with the estimates of Ref . [H] . The proton component 
does not contribute to the Schiff moment of ^^^Hg since the valence nucleon 
is a neutron. However, due to neutron-proton residual strong interaction some 
neutron component is produced by the core polarization. This component is 
shown in Fig. 2. Qualitatively, the behavior of the neutron component inside 
a nucleus resembles the behavior of the bare operator in Fig.l. However, the 
magnitude of the neutron component is smaller and the behavior outside of the 
nucleus is completely different. 



4 Results for ^^^Hg 

In Table 1 we show the contributions of the neutron component discussed above 
to the Schiff moment of ^^^Hg . The three columns correspond to three isospin 
channels. The first row is the contribution from SRdir produced by SUdir given 
by Eq. ITJ. The second row is the contribution from SR^xch produced by SUexch 
given by Eq. 0, and the third row is the sum of them. 

In the direct contributions the isospin channel T = 1 dominates, as it was 
mentioned above. The exchange contributions are small, in general. For the 
isospin channels T = 0, 2 they are comparable to the direct contributions just 
because the Hatter's are suppressed by the factor {N — Z)/A. 

The contributions from 6S are more significant. In Table 2 we list the con- 
tributions from the direct and the exchange parts of the weak interaction W. 
In the first two rows the unrenormalized contributions SSq are shown. Again, 
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Figure 1: The proton component 
of the SchifF moment. Full line is 
the renormalized operator after sub- 
traction of the spurious component. 
Dashed line is the bare operator Eq. 



Figure 2: The neutron component of 
the renormalized Schiff moment. 



Table 1: Contributions of the neutron component to the Schiff moment of ^^^Hg 
in [e fm^]. 





990 


991 


992 


direct 


0.0038 


-0.024 


-0.0076 


exchange 


0.0032 


-0.002 


-0.0015 


total 


0.0070 


-0.026 


-0.0091 



the exchange contributions are considerably smaller than the direct ones. The 
contributions from the channels T = 1, and T = 2 are comparable here although 
T = 1 contribution is still larger. Next two rows show the renormalized contri- 
butions SS. The effect of the core polarization is significant here. The induced 
moment SS includes spin dependent operators, therefore, the spin-spin part of 
the residual interaction Eq. (|16|) is responsible for their renormalization. The 
spin-spin interaction is repulsive with the constants gs — 0.63, and g'^ — 1.01. 
The repulsion results in decrease of absolute values of the renormalized contri- 
butions. Finally, in the last row we show the contribution of the non-local term 
SSnl Eq. I|15|) . It has the exchange origin as well and its contribution is really 
insignificant. All the contributions can be summarized as follows 

S = -0.00045^0 - 0.055^51 -I- 0.009^52 [efm^]- (22) 

The obtained value for the Schiff moment in Eq. (|22() cannot be compared 
directly with the previous calculations with the contact interaction Eq. 
The reason is in different definition of the dimensionless constants ggi in Eq. 
and r\ah in Eq. Q . To perform the comparison we redefine the constants gi 



(23) 



Table 2: Induced contributions to the Schiff moment of ^^^Hg in [e fm"^]. 





990 


991 


992 


direct SSq 


-0.0302 


-0.0631 


0.0604 


exchange SSq 


-0.0007 


-0.0012 


-0.0007 


direct SS 


-0.0086 


-0.0285 


0.0172 


exchange SS 


-0.0002 


-0.0008 


-0.0003 


non-local (5Sjvl 


0.0014 


-0.00004 


0.0013 



With this factor the integration over space of the Yukawa function gives 1, 
exactly as the integration of a S{r). Introducing this factor we obtain 

S = (-O.Olggo - 0.86351 + O.Wggs) x lO^^ [efm^]. (24) 

This value should be compared with = — 1.4 x 10~^rinp from Ref. jJH], and 
S « —1.6 X 10^*77„p from Ref. 6 . Remembering that rjnp ~ g{go + gi — 2^2) 
we conclude that the difference between our result and previous calculations is 
significant for T = 0, and T — 2 channels. Our values are smaller in absolute 
value. In order to trace the origin of this difference we repeated our calculations 
using the contact interaction and omitting the core polarization. The contact 
interaction was obtained by replacing the Yukawa function in Eq. I^J by the 
delta-function. The result is 

S = -0.96 X 10-^gigo + ffi - 2^2) [efm^]. (25) 

If we omit completely the effect of the core polarization in calculations with the 
finite range interaction Eq. 10} then the only nonzero contribution in the first 
row of the Table 2 becomes 

^ = -0.0865(30 + 51-252) [efm^], 

that corresponds to 

S = -1.35 X 10~^5(.go 4- gi - 2.92) [efm^]. (26) 

Comparing Eq. H25|l and Eq. H26|l we conclude that the effect of finite weak 
interaction range is not very significant. The main effect bringing the value of 
the Schiff moment from that in Eq. (|26|l to the value in Eq. H24() comes from 
the core polarization. 

The last remark concerns the pairing effects. The nucleus ^^^Hg has 7 neu- 
tron holes in the unfilled shell. One of them is fixed in pi /2 state and 6 others 
should be distributed among the other states in the shell. From mass differences 
we found A„ — 0.69 MeV. This value of the pairing gap is typical for developed 
pairing. For dipole transitions the transition energy is large compared to A and 
the pairing effects are small. They were omitted in Eq. H17(l . For the induced 
moment the situation is different. There, the transitions with A J = 0, 1 and 
AL — 0, 2 are responsible for the core polarization. Such transitions exist inside 
the last unfilled neutron shell and, due to Pauli blocking, they are sensitive to 
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the details of the shell occupation. For a T-odd operator the effects of pairing in 
the core polarization can be considered in the way used in Ref. jllSj . We found, 
that the main effect of pairing is in fixing the occupation numbers in the upper 
unfilled neutron shell. As soon as we keep the occupation numbers fixed the 
results for SS are changed within few percents when we put A„ — 0. 

In summary, we calculated the Schiff moment of ^^^Hg nucleus using finite 
range weak interaction and considering the core polarization effects. The effects 
of the finite interaction range are not very significant for the Mercury nucleus. 
They do not change the order of magnitude of the Schiff moment calculated 
with the contact interaction. The effects of the core polarization are two-fold. 
First, they renormalize the bare operator of the Schiff moment producing a small 
neutron component due to n-p residual interaction. Second, they produce an 
induced Schiff moment due to P- and T- violation component in the intermediate 
single particle states. The induced moment is proportional to both the strong 
residual and the weak interactions. The effects of the core polarization for the 
Mercury nucleus are large and they have to be accounted in calculations of P- 
and T-violating effects in the Mercury nucleus. 
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